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Abstract 

TeV scale gravity scenario predicts that the black hole production dominates over all 
other interactions above the scale and that the Large Hadron Collider will be a black 
hole factory. Such higher dimensional black holes mainly decay into the standard model 
fields via the Hawking radiation whose spectrum can be computed from the greybody 
factor. Here we complete the series of our work by showing the greybody factors and the 
resultant spectra for the brane localized spinor and vector field emissions for arbitrary 
frequencies. Combining these results with the previous works, we determine the complete 
radiation spectra and the subsequent time evolution of the black hole. We find that, for 
a typical event, well more than half a black hole mass is emitted when the hole is still 
highly rotating, confirming our previous claim that it is important to take into account the 
angular momentum of black holes. 
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1 Introduction 



Black hole has been playing crucial roles for decades in the yet unaccomplished theo- 
retical development to reconcile gravitational interactions with quantum description of 
nature P 13 El IH ISl El El- From this point of view, it is no doubt desirable to have a 
direct experimental test of the Hawking radiation through which black hole is expected 
to radiate off quanta almost thermally. An astrophysical black hole is too big to have a 
sizable Hawking temperature for this purpose. It has long been known that a collision 
of two particles above the Planck energy scale inevitably leads to a production of black 
holes [8j, which can have large Hawking temperature. 

Until recently, the Planck scale was always assumed to be unaccessibly high for human 
beings even in theories with extra compact dimensions such as string theory, under the 
assumption that the compactification radius is of the order of the Planck length. Without 
prejudice, there is no reason that forces one to keep this assumption. Recently it has been 
pointed out that the large extra dimensions can reduce the higher dimensional Planck 
scale down to TeV scale to solve the so-called hierarchy problem without contradicting any 
experimental data if the standard model interactions are confined to a 3-brane [HI El E] • 
It has also been pointed out that a TeV scale gravity is realized if an extra dimension is 
compactified with the warped metric on AdS^ [12]. 

The TeV scale gravity opens up a possibility of producing black holes and observing 
their decay products directly at next generation collider. First it was considered that 
such a black hole will mainly decay into the bulk graviton modes |131 ITl] . Later it has 
been noticed that the decay channels into observable brane-localized standard model fields 
dominate over the bulk graviton emission ISj.When the black hole is highly rotating, the 
bulk graviton emission is expected to be greatly enhanced [IHl HZl CHI UH] • Currently the 
bulk graviton equation is obtained only for the non-rotating J = hole [201 1^ 122] and 
we cannot conclude how much the graviton emission will be enhanced. There is also an 
argument that the higher dimensional graviton has a larger number of modes than four 
dimensions and that the graviton emission would be enhanced due to this fact PHI EI- The 
actual number of modes that have smaller masses than the Hawking temperature, which 
is always smaller than the higher dimensional Planck scale, will depend on the details 
of the moduli stabilization. Typically one has much fewer light modes after the moduli 
fixing. The possibility of black hole production at collider and the detection of its Hawking 
radiation is studied from general perspective in [221 and the experimental signature at the 
CERN Large Hadron Collider (LHC) is studied in |21] in the approximation that the black 
hole has vanishing anglar momentum J = and that the Hawking radiation spectrum is 
described with the high frequency limit ^ 1 , where r/j and u are the black hole horizon 
radius and the energy of the emitted particle, respectively. 

In [25], we have emphasized the importance to take into account the angular momentum 
of the black hole and have first estimated the differential production cross section for a 
given J, whose integral over J qualitatively agrees with the numerical results of the classical 
gravitational collision of two massless particles [2^1 12Z|; namely, the fact that the total 
cross section increases from the Schwarzschild approximation more and more for higher 
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dimensions is well described.^ Further, utilizing the Newman- Penrose formalism jHOj we 
have obtained the field equations of the brane- localized scalar (s = 0), spinor (s = 1/2), 
and vector [s = 1) fields in a separable form for the D = 4 + n dimensional black hole 
background with a single angular momentum J > [SI]. Later these equations are also 
derived for a non-rotating J = black hole in Refs. jS21 EBj- 

Once shown to be separable, its radial part can be solved to obtain the greybody 
factors that completely determine the Hawking radiation spectra without relying on the 
high frequency limit ^1. A thermal radiation spectrum at the horizon is modified 
when the radiation passes through the curved geometry toward an observer located at the 
spatial infinity. The greybody factor is introduced to take this correction into account for 
each angular mode (/, m) and frequency u of the emitted radiation. It is defined to be the 
absorption rate by the hole for a steady in-falling fiux at the infinity 



Im 



Im Im 



(1) 



under the condition that the radial wave is purely ingoing at the horizon. This definition 
physically means, in a time-reversed sense, that F is identical to the proportion that passes 
through the gravitational potential from the horizon toward the infinity without being 
refiected back. 

In early 1970's, Teukolsky showed that the general wave equation in the four dimen- 
sional Kerr black hole background is separable into the angular and radial parts and devel- 
oped analytic and numerical methods to solve the equation ESI 110] . Page then calcu- 
lated the Hawking radiation spectra with the greybody factors taken into account jSZI and 
shown the whole time evolution of the four dimensional Kerr black hole (SHI. The brane 
field equations that we have obtained in |2S] is a generalization of the Teukolsky equation. 
In [22], we have further obtained the analytic formulae of the greybody factors of the s = 0, 
1/2 and 1 brane fields in the low frequency limit rhuj <^ 1 for D = 5 dimensional J > 
black hole. 

Next goal is to obtain the greybody factors without relying on the low frequcy expan- 
sions. The numerical results with non-rotating J = black hole are shown by Harris and 
Kanti in The greybody factors of s = brane field for general J > hole are pre- 
sented by ourselves [1211101 and by Harris and Kanti [Hlll^l. Recently there also appeared 
a paper that treats the s = 1 modes for J > black hole in a restricted range < 4 [IH] . 

In this paper, we complete the series of our work by presenting the greybody factors 
for s = 1/2 and s = 1 brane fields for the general J > hole and the resultant Hawking 
spectra. Now that all the brane field spectra are completed, we are finally able to show 
the whole time evolution of the black hole to confirm our previous claim that the spin- 
down phase, where the approximation J = is not valid, is hardly negligible: we find that 
typically much more than a half of the black hole mass is emitted in the spin-down phase. 



similar formula has appeared in '28', which predicts that the cross section would be reduced by taking 
into account the angular momentum, contrary to our result. See also 1291 for another earlier attempt. 
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The organization of the paper is as follows. In section 2, we go beyond our previous 
treatment to obtain the brane field equations in a separable form for the lower spin com- 
ponents too and recast the results into a rather simple formulae (j5T| and (j52|l . We then 
show the resultant forms of the asymptotic solutions and the greybody factors. In section 

3, we explain our numerical methods to evade the contamination of the outgoing wave 
near the horizon when we impose the purely ingoing boundary condition there. In section 

4, we show our results of the greybody factors and the Hawking spectra for the brane- 
localized spinor and vector fields. In section 5, we combine all the results of scalar [IHI and 
spinor/vector obtained in the previous section to determine the whole time evolution of 
the black hole. The time evolution for the Randall- Sundrum (D=5) and the ADD (D=10) 
black hole are presented. 



2 Brane field equations 

The stationary rotating black holes would be described by the Myers- Perry solution 'WP^ as 
indicated in |^ |^ H^l El HH E]. The induced metric on the three-brane in the (4 + n)- 
dimensional Myers-Perry solution with a single nonzero angular momentum is given by 

A - sin^ ^ , , 2a(r'^ + - A) sin^ ^ , , 
g = dr H ^ dtdip 

Zj Zj 



sin' Mifi' - —dr' - T.d^\ (2) 



S " A 



where 



S = + cos^ ^9, A = + _ ^^i-". (3) 

The parameters /x and a are equivalent to the total mass M and the angular momentum J 

^ ^ (2 + n)A2+„/i J ^ 

evaluated at the spatial infinity of the (4 + n)-dimensional space-time, respectively, where 
A2+n = 27r(3+")/Vr((3 + n)/2) is the area of a unit (2 + n)-sphere and G^+n is the (4 + n)- 
dimensional Newton constant of gravitation. 

In Ref. |2^ it has been shown in terms of the Newman-Penrose formalism that the spin 
s = 0, 1/2 and 1 equations are of separable type in the background space-time with the 
metric (j2I). Here we give the basic equations in more compactified notation. 

The null tetrad {n, n', m, m'} on the brane is 

S 
A' 



n = dt — asin^Wu, — -rdr, 



n' = -^(dt-a sin^ dd^) + ^dr, 

zsint? r 2 2\ n 1 r — ia COS 

v^(r + mcos^9) ^ ^ ^ V2 

m! = m. (5) 
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The spin-coefficients j^O] in this tetrad system become 

ia sin 1 pcot'd 

P= ■■ 0^ P 



a/2S ' r -ia cosd' 2a/2 ' 

mp^sin'i? , P^pA , , PpA^, 



P = e = p 



72 ' 

= r' + /3, K = = = a' = e = 0. (6) 

The null tetrad basis vectors define the differential operators 

D', m, m') = (V„, Vn', V^, V„0- 

In terms of these differential operators and the spin-coefficients, the field equations for 
s = 1/2 fields {ipo,ipi) (Weyl equation) are given by 

D^i-5Vo = (/3'-r')^o + (p-e)^i, (7) 
5^1 -DVo = (e'-p')V^o + (r-/?)V^i. (8) 

On the other hand, the field equations for s = 1 fields (0o;0i;02) (Maxwell equation) are 
given by 

D(j)i-5'(f)o = {2P' - T')(j)o + 2p(j)i - (9) 
D(f)2-5'(l)i = - 2r>i + (p - 2e)02, (10) 

= (p' - 2e')0o - 2r0i + a02, (11) 
D'(l)i-6(j)2 = + 2p>i + (2/3 - r)02. (12) 

The all field variables are assumed to have time and angular dependence 4>a, oc e~*'^*+*™'^. 
Here it is useful to introduce the differential operators 

= d^ + Q + Ncot^, = ~Q + Ncot^ (A^ = 0, ±1/2, ±1) (13) 

where 

K{r) = uj{r'^ + a^) — ma, Q^i!)) = —coasmi!} + mcscd (14) 

have been defined. 

2.1 s = 1/2 field equations 

The s = 1/2 field equations ((Tj) and (jH)) can be written as 

^Air?_i = -^im (15) 

2 2 2 

#A5r/i = ^iVi (16) 



where 

7]_i = V2A-^p-i^,, r]. = iPo (17) 
have been defined. From ^Sf^x Eq. (fT3j) + i^A^x Eq. (fTBj) we obtain 

2 

^A3#A^i = (18) 

2 2 2 '2 

By putting rji = i?i (r)S'i ('i9)e~*'^*"'"*'"'^, we have separated equations 

^A5#A5i?i = Aii?i, (19) 

222 ^ ' 

^l^iSi=-\iSi (20) 

22222 

where Ai is the separation constant. In a similar manner, ^tA2X (Uni) - ^1 X (Uni) and 
the substitution ri_i = _R__i (r)S'_i ('^)e~*'^*"^*™'^ give the separated equations 

^^A^^A^^R_i = XiR_i, (21) 

222 ^ ' 

^1^1 S_i = -XiS^i. (22) 
252 22 

The separation constant must again be A 1 , which can be seen from exphcit forms of angular 

2 

equations and (j221)- Hereafter, the angular function 5'^(?9) is normalized such that 

\SaW sin = ^. (23) 
Furthermore, from Eqs. (fT3j) and (fT^ . we obtain relationships 

S_i = ^^iSi (24) 

' /a7 ^ ^ 



Si = ^L^^5„i, (25) 

' /a7 ^ ^ 



and 

^A^R_i = ^fx^Ri (26) 

2 y 2 2 

#A^i?i = .[\Ir_i (27) 
2 V 2 2 

where the constant of proportionality can be determined by noting that 

fifi^iVsin^? = \iSlsiYi^ + di)(Si^iSi\ (28) 

\22/ 22 \222/ 

(^_j5_i)%in79 = AiS^i sin^9 + 9^ sin79^2^_i) (29) 

hold. 
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2.2 s = 1 field equations 

The s = 1 field equations (j9|)- (fT2|) become 



-1; 



where 



(30) 
(31) 
(32) 
(33) 



(34) 



have been defined. From Eq. (p?n|) — Eq. (jHSI), we obtain 

(^-p)(^^ + p)A + (^o^ — ipa sin -t?) (=2i + ipa sin -i?) 

= (^^# A + ^J-^i - y^) ^1 = 0- 

Putting rji = _Ri(r)S'i('i9)e~*'^*"'"*'"''', we obtain separated equations 

(^#A + 2icur) Ri = XiRi 
[^^^i - 2ujacos^^ Si = -XiSi, 



(35) 

(36) 
(37) 



where Ai is the separation constant. In a similar manner, from Eq. ()31|) — =SfoX 

Eq. we obtain 



(#-p)(^ + p)A + (^o — ipa sin ipa sin-i?) 

f #^A + ^0^1 



^-1 



(3J 



V ' P 

Putting ?7_i = i?_i(r)S'_i('^)e~*'^*+*™''^, we obtain another set of separated equations 

(#^A - 2iujr) = Aii?„i (39) 

(^0=^1^ + 2cjacos?9) = -Ai5_i (40) 

The constant of separation Ai is again common for both sets of equations. 

The relationship between rji and ?7_i is seen as following. From =^oX Eq. — 
Eq. (jnU, we obtain 



From Eq. (jSl - =^o^x Eq. dSl, we obtain 



(41) 
(42) 
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These imply that 



^Q^\S\ — BS—i, 
^^Ai?_i = BRi, 

##Ai?i = BR_i, 
where S is a constant. The constant B is determined by 



— Xl — 4u!a{u!a — m). 



We choose 



B — \J \\ — 'iuja{uja — m). 



2.3 Master equations for brane fields 

The angular and radial equations can be recasted into a neat form 



- AQ + (1 - |s|)coti? j Idi) + + |s|coti? 



sin 7? 



+ A 



s K 



5, = 0, 

s 



+ {2\s\-l)i—K,r-\s\ 

s 



Rs = 0, 



More explicit forms are given by 



sin -d d-d \ d-d 



+ 



{s — Loa cos '&y — {s cot '& + mcsc'&y 



— uja{ijja — 2m) — \s\{\s\ + 1) + \i 



dr V dr 



SsW = 0, 



+ 



- isKA 



RJr) = 0. 
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These expressions are valid also for the massless scalar field by setting s = 0. 
Besides we have relationships between different radial components 



a^ + i^ 1 A^i?i = ./Ali?,! (53) 

A / 2 V 2 2 

for spinor fields and 

(^dr + (^dr + ARi = ^J Xj - 4uja{ua - m)R_i (54) 
for vector fields. 



2.4 Asymptotic solutions 

Asymptotic far field solution to the equation (jK^ is 



R.^ F^)^^ 1-^ 



e ( , .A|,| + s 



' s 



where we have introduced the new radial coordinate 

^ (56) 

and the dimensionless parameters ^ = r/r/^, = Vi^jrh, oj^ = uvh and a* = ajrh have been 
defined. 

Then, Eqs. f)53p and ()54|) give the relationship between coefficients 

2 

A/Af — 4cj^,a*(co'*a* — m) 

2.5 Greybody factors 

The greybody factor is given by the absorption coefficient for the wave equations. By 
virtue of the Eqs. (j37j) and (|^. the absorption coefficient can be calculated solely by 
solving single radial equation. 

The number fiux vector for s = 1/2 is 

f = k (n''?/^i7/^i + n'^ipotpo - mf'-ilJitlJo - m'^tfjoilji) (59) 
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where /c is a constant. The radial component become asymptotically 



2 I |2 
- hi 



4(i^-riv-ir-i^fTi^ir) (r-+oo). (60) 



Therefore, the absorption coefficient is given by 



7\r(out) „ I v(°"*) 1^ 

ri = i-^ = i-^^W, (61) 



2 



4 2 



= 1 




y-^(0Ut) 2 




[\l — 4ci;*a^,(c<j^,a^, 


-m)] y^('°) ^ 



for spinor fields, where iv[™^°"*'* denotes the total number flux. The last expression is 

2 

determined by solving only s = 1/2 equation. 

In a similar manner, we have the absorption coefficient 



ri = i-^ = i- ^.2 ^M,.:...,2 > (62) 

for vector fields. 

3 Numerical methods to obtain greybody factors 

We explain our numerical methods. In this section, we take the unit 

Th = 1 (63) 

and always consider the case s > unless otherwise stated. First, we switch from the 
Boyer-Lindquist frame Q to the ingoing Kerr-Newman frame by: 

dv =dt + ^ dr (64) 
dip =dip + (do) 

The radial wave equation now becomes 

d^R dR ^ „ , , 

_ + ,,_ + Cfl = 0. (66) 

with 

is - 1)A' + 2iK ^ 2iujr(2s - 1) - A 

C = X (67) 



A ' ^ A 
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For the angular eigenvalue A, we use the expansion of the form 

A = ^C„(a^r (68) 

n 

presented in jSHj. The coefficients C„ are given for each angular mode and dumps 

quite fast, like a multiple of inversed factorials, as / increases. We take up to n = 6 terms. 
The accuracy of the approximation can be checked by the relative size of the last n = 6 
term to the sum when au is within the relevant region to the peak of the power spectrum, a 
posteriori. We have checked that the ratio is at most of order .1% for each relevant region. 

The asymptotic solutions in the near horizon (NH) and far field (FF) regions are ob- 
tained as (2^ 

R'''' ~ IVin + Wonte^'^'-A^ i?^^ ~ Y^y^-' + Foute'^'^^V^, (69) 

where r* is the tortoise coordinate defined by dr^,/ dr = (r^ + a^)/A and r* — > r for r ^ oo. 
Recall that A ^ for r 1. This coordinate change is not essential for the following 
analysis but convenient because the near horizon ingoing solution does not contain the 
tortoise coordinate r^,. 



3.1 Removing outgoing contamination at near- horizon 

For the case of scalar (s = 0) the calculation is simply to put the purely ingoing 
boundary condition Wi^ = at a point tq = 1 + e close enough to the horizon e <C 1 and to 
solve the second order ordinary differential equation ()66|1 numerically toward the far field 
region. In this region we can easily read off the coefficients l^n,^ut by the fit, whose 
ratio directly leads to the greybody factor F = 1 — jl^ut/^n In- 
putting the purely ingoing boundary condition Wi^ = at the near horizon r = rg is 
always polluted by a tiny outgoing wave, numerically. This itself is the case for the scalar 
too. The difficulty in the spinor (s = 1/2) and vector (s = 1) case is that the amplitude of 
the outgoing wave grows with respect to that of the ingoing one as we go apart from tq. 
To evade this problem, we first expand the near horizon solution 

RfJ^ = 1 + ai(r - 1) + a2(r - 1)' + ■ • • , 

Rout = e'^'^'-ir - 1)^ (1 + 6i(r - 1) + . . . ) , (70) 

where the coefficients aj, bi are straightforwardly obtained by substituting the expansion 
of A,r],( and (f7n|) into (jU^ . which serves linear equations for the coefficients. The point 
to remove the outgoing contamination is the following subtraction 

R = R-{l + ai{r-l)), (71) 

introduced in as Bardeen's prescription. Then R satisfies the equation 

CR = g, (72) 
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where 



£ = (f jdr"^ + rjd/dr + (, 

g = -C (1 + ai(r - rn)) = -rj^i - C (1 + - th)) ■ 



(73) 



Now (ff^ can be safely solved toward the far field region without the growing outgoing 
contamination. 



3.2 Matching at far field 

We expand the far field solutions for spinor (s = 1/2) 



FF 



FF 



/ / / 

1 + — + 4 + 4 + - 



d{ d{ 



- 1 + ^ + ^ 



(74) 



and for vector (s = 1) 



R. 



FF 



«f «u r'" 

rfl + ^ + % + % + - 



R. 



FF 



(75) 



Again the coefficients cf'^, d{'^ can be straightforwardly obtained. The asymptotic form of 
the i?, is now 



R 



1/2 



(l^n + 1 - ai) + air + l^n^ + ^'out^— , 



Ri = (lin - ai) r + (l^ncl - 1 + ai) + l^n- + >;ut- 

r 



2iajr* 



(76) 
(77) 



By the fit to the obtained solution R in the far field region, we can determine Y^^, Fout 
since the bigger terms than that contains yout are all fixed. The smallness of the oscillating 
term containing Yont can be easily overcome by keeping sufficient digits in the numerics 
and by taking sufficiently dense sample points in the far field region for the fit- 
Finally the greybody factors can be obtained as 



s=l/2 



1 



2ui 

"RT 



Yin 
Yout 



Yi„ 



(78) 
(79) 
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4 Greybody factors and Hawking radiation spectra 



In j2Sl we have argued that the black hole production cross section for a center of mass 
energy s is well approximated by 

^( \ - SttJ/s for J < Jmax, /o^N 

dJ^'^ forJ>J^,., ^""^^ 

except for the region where J is very close to Jmax, and consequently that the black 
hole tends to be produced with large angular momentum. The rotation parameter a 
corresponding to J^ax is |2S1 

71 

flmax = 2~^^ ^^^^ 

for D = 4 + n dimensional black hole. 

Once produced, we assume that the decay process of the black hole is governed by the 
Hawking radiation into the brane-localized standard model fields See |25 for a review on the 
estimation of the amount of energy radiated at the black hole formation process (balding 
phase) rather than by the Hawking radiation and see also [221 for the recent progress, in 
the spirit that the quantum gravitational correction will be read off as a deviation from 
this precise prediction in the black hole picture. 

With this assumption in mind, the number of spin s particles emitted into an (Z,m) 
angular mode of the spheroidal harmonics is given, per degree of freedom per unit time 
per energy interval [u, u + duj], by 



2s' 



where T = [(n + 1) + (n — l)a^]/47rr/t(l + al) and Q = a*/(l + al)rh are the Hawking 
temperature and angular velocity of the black hole, respectively. The corresponding power 
and angular spectra are obtained by multiplying the number spectrum (jH^ by u and m, 
respectively. 

The time evolution of the black hole with mass M and angular momentum J is then 
governed by 

where Qs is the number of massles degrees of freedom at temperature T, namely, the number 
of degrees of freedom whose masses are smaller than T with spin s. 



4.1 Greybody factors for spinor and vector fields 

We present the greybody factors obtained by the procedure above. Hereafter we limit 
ourselves to the case D < 11 motivated by the fact that D = 11 is the highest possible 
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dimension if we exclude the s > 2 component fields from the supergravity multiplet, though 
in principle we can consider a larger dimensional case as well |53j . 

In Fig. 1, we plot the greybody factors for the brane- localized spinor field for the non- 
rotating (a* = 0) and highly-rotating (a* = 1.5) black holes. Note that a^: = 1.5 is the 
highest possible rotation parameter for a D = 5 black hole when the black hole is produced 
by a collision of two particles, see Eq. (|HT|). 

The Figure 2 shows the corresponding plots for the vector field (s = 1). We observe that 
the greybody factor becomes negative when uj < mQ, rendering the number spectrum (jH^ 
always positive. In other words, when the black hole is highly rotating, an incoming steady 
energy flux with u < mQ is scattered back by the hole with an increased amplitude, see 
Eq. . This is called the superradiance [53] . Though we cannot observe the superradiance 
itself for the TeV black hole since it is decaying so fast, we show the rate of the energy 
amplification in Fig. 3 for its own interest. 

4.2 Power and angular spectra for spinor 

Hereafter in this section, we present the dimensionless power and angular spectra per 
degree of freedom: 

^'^ dudt ~ 27re(^-'"f^)/^ - (-1)2^' dudt ~ 2tt ei"^)/^ - (-1)'^'' ^ ^ 

versus the dimensionless energy of the emitted particle VhU. Note that this implies that 
the frequency cu in the horizontal axis is given in unit of r^^^, which varies for a fixed mass 
M when we vary the angular momentum J. For simpler presentation we show our plot for 
< ThUJ < 8 in this section but we have also calculated all the u regions for a* < 1.5 to 
obtain the total power and angular spectra in the next section, see also Figs. 14. 

In Figs. 4 and 5, we plot them for the brane-localized spinor and vector fields, respec- 
tively, with varying black hole rotation a^, for fixed D = 5 (Randall-Sundrum black hole) 
and D = 11 (ADD black hole). It is clear that the larger the rotation parameter a* is, the 
more enhanced both the power and angular spectra are. 

In Figs. 6 and 7, we plot the same for spinor and vector, respectively, for varying 
dimensions D with fixed black hole rotation a* = (non- rotating) and a* = 1.5 (highly- 
rotating). These figures show that the black hole radiation is greatly enhanced for larger 
dimensions. Note that the peaks coming from e = m modes are distinctive only for the 
Randall-Sundrum black hole, especially in the case of the spinorial emission. 

The substructure behind the total spectrum is shown in Figs. 8,9,10 and in Figs. 11,12,13. 
For each column, the greybody factor, power spectrum, and the angular spectrum are 
shown from above to below. We can see that each region of the greybody factor's rise in 
coincides with each peak. For non-rotating hole, the angular modes m = —I, —l + l,...,l 
give exactly the same greybody factors and the ±m modes cancel each other in each 
given / mode for power and angular spectra. When the black hole is highly rotating, the 
contributions from the I = m modes become dominant. 
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For a highly rotating black hole with large a^, the contributions from / = m modes 
become dominant, which greatly enhance the amplitude of total spectra at high frequency 
region rhUJ > 1. This is sometimes called the superradiant enhancement of the higher spin 
emission, though it is not directly related to the original meaning jHll of the superradiance 
as shown in Fig. 3 To exhibit such a high frequency tails which are taken into account in 
the calculation of the next section, we present Figure 14. 



5 Time evolution of black hole 

In Fig. 15, we show the schematic pictures for the black hole life at various stages: pro- 
duction phase, the balding phase, the spin-down phase, the Schwarzschild phase, and the 
Planck phase The spin-down and Schwarzschild phases are of interest here. 

We show the black hole time evolution as it emitts brane localized particles. We calcu- 
late the rates at which energy and angular momentum are radiated into the brane-localized 
standard model fields and the evolution of the mass and the angular momentum. A partic- 
ular quantity of interest is the portion of the energy emitted during the spin-down phase 
and the (almost) Schwarzschild phase. 



5.1 Setup for time evolution 

To follow the time evolution, it is convenient to make quantities invariant under the scaling 
of rs = {1 + aiy^^'^'^^^rh, which does not vary for a fixed M when we vary J. We define 
a scale invariant function 7(05), with respect to the scale invariant rotation parameter 
ttg = a/vg, as follows: 

7 = ,1 ,^/,^ (85) 



d\nM/Mi 
n + 2 fl dJ 



a dM n 



^6) 



where the mass of a hole is measured in the unit of the initial mass Mj. This quantity is 
directly integrated. We calculate the ratio of the final mass Mf to the initial mass Mj by 
integrating Eq. ()85|) with an initial rotation parameter a5(ini) 

-/ = exp / da,^ . (87) 

The amount of energy which is radiated in spin-down phase (0 ~ (final) ^ ^ as(ini)) 
is {Mi — Mf) and then the remaining Mf will be subsequently radiated off in the Schwarz- 
schild phase, where the angular momentum of black hole is vanishing. 

We consider the evolution of the black hole. Since the time roughly scales as r""*"^ in 
(4-|-n) dimensions, it is convenient to introduce scale invariant rates for energy and angular 
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momentum as follows.^ 

,d\nM 



a(a^ = -r"+3 



dt 

Pias) ^ --r'^' (89) 

with these new variables 7(0^) can be written as ■y~^{as) = P/a{as) — (n + 2)/(n + 1). For 
all the standard model particles, 

asM = 93(^3=0 + gf<ys=i/2 + gvas=i, (90) 

PSM = gsPs=0 + gfPs=l/2+ 9vPs=l, (91) 

where (7^ = 4, (7/ = 90 and = 24 are adopted in this section. 

We also introduce dimensionless variables y and z to take angular momentum and mass 
of the hole: 

y = -lna„ (92) 
M 

. - ->"-5i. (93) 

then finally we get the time variation of energy and angular momentum in terms of scale- 
invariant time parameter (r = r~"~'^(ini)t) with initial mass of the hole: 

dz a 



dy /5-«(Sf)' 

-f = (/?-« ^ ^— ) e"+l^ (94) 



dr \n + 1 

After finding the solutions z{y) and r(?/) of the coupled differential equations OH one 
can get |/(r) and z{t\ hence and M/Mi, as a function of time. From these, one can find 
how other quantities evolve, such as the area. 

For our purpose, it is convenient to convert the set of variables (r/j, a^) For conversion 
of unit, the following expressions are useful with = a*/(l + 0^)^/*^"+^^ 

«(a,) = -C\l + al)^^rl^, (95) 
P{as) = -<^'(l + «*)^^^^^' (96) 



where 



TsM-^^ = ——— , (97) 



[n + 2)a 



n+2 



,n + 2^ 1 

(9^ 



^We can understand this by simply looking at the formula —dM/dt ^ AT'^ where the surface area of 
horizon A ^ for brane fields and the temperature of the hole T ^ l/r^ and M ^ 
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5.2 Results for time evolution 



In this section we use the natural unit SirG = 1. We assume that all the standard model 
fields are massless and therefore the effective degrees of freedom are given by Qs = 4, 
(?/ = 90 and = 24. In Fig. 16, we draw the rates at which the energy and angular 
momentum are radiated into the brane localized standard model fields. We have explicitly 
calculated the rates up to the rotation parameter a* < 1.5 for the Randall- Sundrum black 
hole {D = 5) and for the ADD black hole {D = 10). For D = 10 we have extrapolated our 
result using the cubic-curve approximation up to a* < 4 to cover all the possible rotation, 
see Eq. (18111 . The curves are plotted with respect to rather than a^, using the conversion 
rules described in the previous section, namely up to < 0.83 and 2.67 for D = 5 and 10, 
respectively. The D = 10 results are exact without using the cubic-curve extrapolation up 
to as < 1.27. 

In Fig. 16, one can see that the power a and torque (3 are increasing functions of the 
angular momentum a^. The vector emission dominates over the spinor and scalar emission 
for the high rotation parameter, sometimes called the superradiant enhancement of the 
higher spin particle emission, but as rotation becomes slower the fermion channel becomes 
increasingly important. Generally, angular momentum is emitted much faster than energy, 
therefore rapidly rotating black hole spins down to a nearly non-rotating state before its 
mass has been radiated off completely. 

In Fig. 17, we plot the mass of a hole as a function of the rotation for the virtual 
setup where only the scalar (s), spinor (f), or vector (v) field is emitted, respectively, and 
for the realistic case where all the standard model fields are emitted (SM). One can see 
that the larger the particle's spin is, the more effectively the black hole angular momentum 
is carried away. For the most effective case of vector-only (v), the angular momentum is 
carried away so rapidly that more than 30% of the mass still remains after the spin-down 
phase, to be radiated off at the Schwarzschild phase. In contrast, the scalar-only case 
(s) exhibits that the black hole radiates its whole mass before it stops rotation. For the 
realistic case (SM), black hole loses roughly 70% to 80% in Z) = 5 and D = 10, respectively, 
before it stops rotation when starting from the maximum rotation. Note that the scalar 
emission (s) is subdominant comparing to the vector (v) and spinors (f) because of its 
small effective degrees of freedom and small emission rates. 

In Figs. 18 and 19, we plot the time evolution of rotation parameter and mass The unit 
time to is defined by the time duration from the initial state to the state with = 0.01 for 
the virtual cases of spinor-only (f) and vector-only (v). For scalar-only case (s), the mass 
goes to zero before rotation stops, therefore we defined to for scalar to be the duration 
until the whole mass is radiated off. The initial rotation parameter is fixed by = 0.83 
and 2.67 in D = 5 and D = 10, respectively. The mass of the hole goes to zero before the 
rotation parameter goes to zero when only scalar emission is available. We have taken the 
initial radiation parameter to be = 0.83 and 2.67 for D = 5 and 10 that are maximally 
rotations allowed by the initial collision. 

When all the standard model fields are turned on (SM), the evolution is essentially 
determined by the spinor and vector radiation. The figures show that a black hole spins 
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down quickly at the first stage with large rotation parameter and the decrease of rotation 
parameter slows down as angular momentum of the hole is reduced. 

6 Summary and Discussion 

We have generalized our previous result [23] for the greybody factor and Hawking radiation 
spectrum of the higher dimensional Kerr black hole to general dimensions D = 4 + n, 
without relying on the low frequency limit Vhu: <^ 1, not only for the scalar field [s = 0) jlHI 
but also for the spinor [s = 1/2) and vector fields {s = 1). Now we can completely describe 
the evolution of a black hole with any given initial mass and angular momentum by taking 
into account any type of field residing on the brane. 

We have developed the numerical method to solve the radial Teukolsky equation which 
has been generalized to the higher dimension {D = 4 + n) in the first paper [25j of our 
series of publication. There are two points in our numerical methods. First, we have im- 
posed the proper purely-ingoing boundary condition near the horizon without the growing 
contamination of the out-going wave by extracting lower order terms explicitly. Second, 
we have developed the method to fit the in-going and out-going part from the numerically 
integrated wave solution at far field region by explicitly obtaining the next-to-next order ex- 
pansion (or next-to-next-to-next order in vector case) of the solution. With these progress 
in numerical treatment, we can safely integrate the generalized Teukolsky equation up to 
very large r without out-going wave contamination. 

Then we have calculated all the possible modes to completely determine the radiation 
rate of the mass and angular momentum of the hole. Totally 3407 modes are computed 
explicitly, other than the modes which are confirmed to be negligible. A black hole tends 
to lose its angular momentum at the early stage of evolution. However the black hole still 
have a sizable rotating parameter after radiating half of its mass. Typically more than 70% 
or 80% of black hole's mass is lost during the spin down phase. In the case of very fast 
initial rotation, the number could be modified quantitatively by taking the bulk graviton 
emission into account, especially for a larger number of extra dimensions as discussed in 
the Introduction, but the result would remain the case qualitatively.^ 

We have determined the radiation and evolution of the spin-down and Schwarzschild 
phases up to the ambiguities for the initial fastest rotations shown above. The remaining 
hurdle is the evaluation of the balding phase, which is still being disputed due to its non- 
purturbative nature, to extract the quantum gravitational information at the Planck phase 
from the experimental data at LHC 
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Figiin; 3: Siip(!rra,dia.nt. Jimplification (if bramvlocalw^d v<x:tor fiold a.mplitiid<^ s<;att(^rod by 
a, highly rotating (a.^ = 1.5) bbudt hole in D = 5 (h^ft) a.nd 10 (right). 
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Figun; 4: Spinor power (iippnr) and angular (lowor) Hp<x:tra for fix<^d D = 5 [h^ft) a,nd 
D = 11 (right). In <^a.(;li figure, a^^ = 0. 0.3, ... , l.T) from lowcr-h^ft (r<xl) to upper-right 
[ma,g<^nta). 



24 



powec spectirum [D=5J powe: spectirum [D=l 1) 
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(magdnrti). 




Figure 6: Spinor power (upper) find Miignhi.r (lower) Kpectra fetr fixed f) (left.) and 1.,') 
(right). Ill (ifieli figiirf;. D 7. 9 tmd 11 frcmi b(4ow ro tiibov*;. 
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Figiiic 7: Vector power (upper) find fiiignlar (lower) Kprcrra for fix(^d a■^: f) (left.) and !.,'> 
(rig'ht.J. Ill efitii figure. D 7, 9 fmd 11 from below to ribfiv*^ 




Figiu-c :S : Spmor greybody f^u;tor ior [D,a,] ( ^i, (]) , {5,1.5). ( 1 1 J] ) fmd (11^1.5). The 
eoiitribiitions from the aiigulai mfjdf^s with the Hni.me / mod 10 ni.r(^ dra.wii with thc^ stime 
color, 
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l?igiitc !l: Spinof pow('r spcx'tfii for {D.a.,) ('),()), (r),l.r)), (11,0) imd (ll.l.rj). The 
t:("tiitribnt:ions fnmi the, aiigulai mfjdos with the Ha.mr I mod 10 uro dra.wii with t.ht^ same 
eokir, 




Figiu-c 10: Spinor ;>ngiila.r spcctm for (D,a,) (r,J)), (r, a. (11; 0) ;md (11,1.5). The 
eoiitribiirions frrtm the aiigulai mfjdf^s with the Humr / mod 10 ni.f(^ dra.wii with thc^ same 
color, 
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Pigiitc 1 1 : Vector grcybody fa,(;roL for {D,a,] ( 5 Jl ) . ( fj , 1 . fj ) . ( 1 1 , 0) and ( 1 1 , 1 ■ 5) . The 
t:("tiitribnt:ion?^ from the, angular mfjdos with the Ha.mr I mod 10 a.rc? drjxwii with t.ht^ same 
eokir, 
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Figure 12: Vector power spectra for {D.aJ) ('j.d), (5.1.5), (11.0) ;uid (11,1,5). The 
eoiitribiitions frrtm the aiigulai mfjdf^s with the HH-me / mod 10 ni.r(^ drawn with thc^ same 
color, 
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Figure 13: Vector angular specl.ra for (D.a*) — (5,0), (o;1.5), (11.0) and (11,1.5). The 
coiil.ribuLioii& from llic aiigiilar modes wil.h Ihe same I mod 10 arc drawn willi llie same 
color. 




Figure 14: Power spectrum Irom D = 10 dimensional black liole for spinor (left) and vector 
(right). Left: similarly plotted as in Figs. 8, 9. 10. Righl : contributions ol I — rn modes 
are shown. 
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